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Abstract 

This paper considers a multiple-input multiple-output (MIMO) Gaussian wiretap channel model, 
where there exists a transmitter, a legitimate receiver and an eavesdropper, each equipped with multiple 
antennas. In this paper, we first revisit the rank property of the optimal input covariance matrix that 
achieves the secrecy capacity of the multiple antenna MIMO Gaussian wiretap channel under the average 
power constraint. Next, we obtain necessary and sufficient conditions on the MIMO wiretap channel 
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I. Introduction 

The broadcast nature of a wireless medium makes it very susceptible to eavesdropping, where the 
transmitted message is decoded by unintended receiver(s). Recent information-theoretic research on 
secure communication has focused on enhancing security at the physical layer. The wiretap channel, 
first introduced and studied by Wyner [1], is the most basic physical layer model that captures the 
problem of communication security Wyner showed that when an eavesdropper's channel is a degraded 
version of the main channel, the source and destination can achieve a positive secrecy rate, while ensuring 
that the eavesdropper gets zero bits of information. The maximum secrecy rate from the source to the 
destination is defined as the secrecy capacity. The Gaussian wiretap channel, in which the outputs at the 
legitimate receiver and at the eavesdropper are corrupted by additive white Gaussian noise, was studied 
in El. 

Determining the secrecy capacity of a Gaussian wiretap channel is in general a difficult non-convex 
optimization problem, and has been addressed independently in Gl-EI. Oggier and Hassibi [3] and Khisti 
and Wornell H followed an indirect approach using a Sato-like argument and matrix analysis tools. They 
considered the problem of finding the secrecy capacity of the Gaussian MIMO wiretap channel under 
the average total power constraint, and a closed-form expression for the secrecy capacity in the high 
signal-to-noise-ratio (SNR) regime was obtained in [4]. In [5], the rank property of the optimal input 
covariance matrix for the secrecy rate maximization problem is discussed but the authors were unable 
to characterize the solution for the general case. For some special cases of the MIMO wiretap channel, 
where the solution has rank one, the optimal input covariance matrix that achieves the secrecy capacity 
under the average total power constraint was obtained in 

In lUl, Liu and Shamai propose a more information-theoretic approach using the enhancement concept, 
originally presented by Weingarten et al. IflOll . as a tool for the characterization of the MIMO Gaussian 
broadcast channel capacity. Liu and Shamai have shown that an enhanced degraded version of the channel 
attains the same secrecy capacity as does a Gaussian input distribution. From the mathematical solution 
in 1H it was evident that such an enhanced channel exists; however it was not clear how to construct such 
a channel until the work of [9], which provided a closed-form expression for the secrecy capacity under 
a covariance matrix power constraint. While this result is interesting since the expression for the secrecy 
capacity is valid for all SNR scenarios, there still exists no computable secrecy capacity expression for 
the MIMO Gaussian wiretap channel under an average total power constraint. 

In this paper, we first investigate the rank property of the optimal input covariance matrix that achieves 
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the secrecy capacity of the general Gaussian multiple-input multiple-output (MIMO) wiretap channel 
under the average total power constraint, where the number of antennas is arbitrary for both the transmitter 
and the two receivers. Next, we obtain the optimal input covariance matrix for the case that this optimal 
covariance matrix is full-rank. Necessary and sufficient conditions to have a full-rank optimal input 
covariance matrix are characterized as well. 

The rest of this paper is organized as follows. In the next section, we describe the assumed mathematical 
model and revisit the current solution for the wiretap channel under the matrix power constraint. The 
rank property of the optimal input covariance matrix under the average power constraint is investigated 
in Section III, and in Section IV we characterize the conditions under which the input covariance matrix 
that achieves the secrecy capacity of a wiretap channel under the average power constraint is full-rank. 
In Section V, we discuss some interesting facts regarding the optimal solution, and in Section VI we 
present numerical results to illustrate the proposed solutions. Finally, Section VII concludes the paper. 

Notation: Vector-valued random variables are written with non-boldface uppercase letters (e.g., X), 
while the corresponding non-boldface lowercase letter (x) denotes a specific realization of the random 
variable. Scalar variables are written with non-boldface (lowercase or uppercase) letters. The Hermtian 
(i.e., conjugate) transpose is denoted by (.) H , the matrix trace by Tr(.), and I indicates an identity 
matrix. Inequality A ■< B means that A — B is Hermitian positive semi-definite. The Euclidean norm of 
the vector x is written as ||x||. Mutual information between the random variables A and B is denoted 
by I(A;B), E is the expectation operator, and CAf(0, a 2 ) represents the complex circularly symmetric 
Gaussian distribution with zero mean and variance a 2 . 

II. System Model and Prior Works 

We begin with a multiple-antenna wiretap channel with n t transmit antennas and n r and n e receive 
antennas at the legitimate recipient and the eavesdropper, respectively: 

y r = Hx + z r 

(1) 

y e = Gx + z e 

where x is a zero-mean nt x 1 transmitted signal vector, z r G C TlrXl and z e G C n ° xl are additive 
white Gaussian noise vectors at the receiver and eavesdropper, respectively, with i.i.d. entries distributed 
as CJ\f(0, 1). The matrices H G c nrXra * and G G C n<,xn < represent the channels associated with the 
receiver and the eavesdropper, respectively. Similar to other papers considering the perfect secrecy rate 
of the wiretap channel, we assume that the transmitter has perfect channel state information (CSI) for 
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both the legitimate receiver and the eavesdropper. For the Gaussian channel, where Gaussian inputs are 
an optimal choice, the secrecy capacity is given by 



C sec = max[J(X; Y r ) - I(X; Y e )] = maxfl(Q) 
x QhQ 



(2) 



where R(Q) = log |HQH^ + 1| — log |GQG^ + 1|, and Q = £{xx fl } is the input covariance matrix. 

In (5J, the above secret communication problem was analyzed under the matrix power-covariance 
constraint, defined as 



Q < S 



(3) 



where S is a positive semi-definite matrix. An explicit expression for the secrecy capacity under Q was 
obtained via applying the generalized eigenvalue decomposition to the following two positive definite 
matrices 



(S^H^HS^+I , S^G^GSs+I) 
In particular, there exists an invertible generalized eigenvector matrix C such that |[T4l 



SfG^GSt +1 



S^H^HSs +1 



A 



(4) 



(5) 



(6) 



(7) 



where A = diag{Ai, A rtt } is a positive definite diagonal matrix and Ai, A„ t represent the generalized 
eigenvalues. Without loss of generality, we assume the eigenvalues are ordered as 

Ai > ... > X b > 1 > A 6+1 > ... > A nt > 

so that a total of b (0 < b < n t ) are greater than 1. Hence, we can write A as 

Ai 
A 2 

where Ai = diag{Ai, A;,} and A 2 = diag{Afo + i, A„ t }. We can partition C similarly: 

C = [d C 2 ] (8) 

where Ci is the n t xb submatrix representing the generalized eigenvectors corresponding to {Ax, A&} 
and C 2 is the n t x (n t — b) submatrix representing the generalized eigenvectors corresponding to 
{A&+1, A nt }. Using the above notation, the secrecy capacity of the MIMO wiretap channel under 
the matrix power constraint ([3]) can be expressed as (5J, lTT3l Theorem 3]: 
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Corollary 1. Under the matrix power constraint ([3]), the secrecy capacity of the MIMO Gaussian wiretap 
channel is given by 



C sec (S) = ^2 lo S A * = lo § I A i I 

1=1 

where the optimal input covariance matrix Q* s that maximizes ^ and attains ([9]) is given by 

(Cfd)- 1 



(9) 



Q* s = S^C 











(10) 



Remark 1. From §5j and one can easily confirm that if H^H ^ G H G, then for any S y we 
have A X I. In other words, in this case the pencil in ([4]) has no generalized eigenvalue bigger than 1. 
Thus, C sec (S) = for any S h 0. 

In this paper, we consider the secrecy capacity problem in (|2]) under the average power constraint: 



Tr(£{xx H }) = Tl-(Q) < P t . 



(11) 



For this constraint, no computable secrecy capacity expression has been derived to date for the general 
MIMO case. In principle, one would have to find the secrecy capacity through an exhaustive search over 
the set {S : S h Q,Tr(S) < P} flU Lemma 1], El: 



C S ec(Pt) = max C S ec(S) . 

ShO,Tr(S)=P t 



(12) 



where for any given semidefinite S, C sec (S) should be computed as given by ^. 

In the next section, we investigate the rank of the optimal input covariance matrix Q* that attains 
C sec (Pt). Next, in Section IV, we obtain the optimal Q* under the average power constraint for the case 
that Q* is full-rank. 



III. Rank Property of the Optimal Solution under an Average Power Constraint 
First, we note that the problem under Tr(Q) < P t is equivalent to that under Tr(Q) = P t (3j , |[5lF 



Also note that in ( 12 1, this implies that we have Tr(S) = P t instead of Tr(S) < P t . 



We are interested in finding the optimal S which maximizes the problem (12i. Let us assume that 



we have found the optimal S. Consequently, from (10 1, the optimal input covariance matrix that attains 



For this statement, and also for the following results in the paper, we exclude the special case H H H -< G H G for which 
the Csec is trivially for any S >z 0, and consequently for any Pt, as pointed out in RemarkJTj 
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C sec (P t ) is given by 



Q* = S2C 



(Cfd 








c"s 



(13) 



where C and Ci have respectively the same definitions as those of C and Ci, given by (|5j)-(|8j, but here 
for the pencil (S^H^HS^ + 1 , S^G^GS^ + 1). Note that Q* can be rewritten as 

(CfCi)- 1 


§s Ci(CfCi) -1 Cf ss 



Q* = S = 





"cf " 







S5 Pg Si 



where Pg = Ci(C^Ci) C± is the projection matrix onto the space of Ci. Moreover, let P 
I — Pg be the projection onto the space orthogonal to Ci. We have 



(14) 



Tr(Q*) = Tr(S^ Pg S*) 
= Tr(SP 6i ) 
= Tr(SP e P 6i ) 
= Tr(P ei SP 6i ) 



(15) 
(16) 
(17) 



where ( 15 ) comes from the fact that Tr(AB) = Tr(BA), and ( 16 1 is because Pg = Pg Pg Similarly 



we have 



Tr(S) = Tr((Pg i +pi i )S(Pg i +Pi i ; 
= Tr(Pg i SPg i ) + Tr(pi i SPi i ) 
= Tr(Q|) + TV(Pi i SPi i ) 



(18) 
(19) 



where ( |T9| > results from (jTT L 

Lemma 1. For the optimal S, we have span{Ci} = span{S}. 



Proof: : The proof is obtained using (19 1, and by noting that for the optimal S we must have 



Tr(S) = Tr(Q*) = P t . This means that we must have Tr(Pi SPi ) = 0, or equivalently Pi S = 0, 

Ci Ci Ci 

which completes the proof. ■ 
Using Lemma [T] in ( p"4| ) we have 



Q* 



(20) 
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The following lemma reveals another property of the optimal input covariance matrix under the average 
power constraint. 

Lemma 2. For the optimal S, i.e. Q*, the pencil (S^H^HS^ +1 , SaG^GS^ +1) has no generalized 
eigenvalue less than one: 

f Ai 



S^H^HSs +1 



S^G^GS^ +1 



I 

I Q 
I 



(21) 



(22) 



where A2 = I corresponds to the generalized eigenvalues equal to one. 

Proof: We note that any vector which lies in the null space of S can be a generalized eigenvector 
of the pencil (S^H^HSs + I , SaG^GS^ + I), with a generalized eigenvalue equal to 1. Such 
vectors span the null space of S, i.e., span{C2} = spanjS} -1 ". On the other hand, from Lemma [l] 
span{Ci} = span{S}. Thus for the optimal S, all generalized eigenvectors of the pencil (S^H^HSa + 
I , SsG^GSs + 1) correspond to generalized eigenvalues either bigger than or equal to 1. ■ 
Let b denote number of generalized eigenvalues of the pencil (S^H^HSs + I , S^G^GS^ + I) 
that are strictly bigger than 1, where again S = Q* represents the optimal input covariance matrix that 
attains the secrecy capacity under the average power constraint given by ( fTT) . From Lemma [T] we have 
rank(Q*) = rank(Ci) = b. 



Theorem 1. For the optimal Q* we have 

rank(Q*) < m 

where m is the number of positive eigenvalues of the matrix H^H — G^G. 



(23) 



Proof: Subtracting ( [22] ) from ( [2T[ ), a straightforward computation yields 

Ai-I 



§5 [H H H — G^G] §2 = C~ H 







C" 1 >- 0. 



(24) 



From (24 1, we note that §2 [H H H - G^G] §2^0, from which it follows that rank(Q*) < 



in. 



Remark 2. From Theorem [T] one can easily confirm that the optimal Q* can be full rank only in the case 
that m = rit, H H H >- G H G. For all other scenarios, the optimal Q* will be low rank. The authors 
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of 0, (51 use the Karush-Kuhn-Tucker (KKT) conditions on problem (|2]) to make a similar statement, 
but they do not show what the rank of the optimal Q* will be. 

The following lemma will be used for the computations in the next section. 

Lemma 3. For the case of H^H y- G H G, for any n t x n t matrix S >z 0, all the generalized eigenvalues 
of the pencil (S» H^HS^ + 1 , G^GSs + 1) are strictly bigger than 1, i.e. Ay I, iff Sis full rank, 
i.e. S >- 0. 



Proof: The claim is easily proved by considering the rank of both sides of (24). 



IV. Characterization of the Optimal Full-Rank Solution 

In this section, we characterize the secrecy capacity under the average power constraint for a particular 
class of MIMO Gaussian wiretap channel where the optimal solution Q* is full rank. While necessary 
conditions for a full-rank Q* were characterized in the previous section, here we derive sufficient 
conditions as well. 

We begin by rewriting problem ([2]) here: 

C sec (Pt)= max log|HQH H + I| -log|GQG H + I| . (25) 

QbQ, Tr(Q)=P t 

The Lagrangian associated with this problem is given by 

£ = log |HQH H + I| - log |GQG^ + I| - /x(Tr(Q) - P t ) + Tr(MQ) (26) 

where \i > and M ^ are the Lagrange multipliers. The optimal Q* must satisfy the following KKT 
conditions: 

H^(HQ*H^ + I) -1 !! - G H (GQ*G H + IJ^G = /xl - M (27) 

//(THQ*) - Pt) = (28) 

Q*M = MQ* = . (29) 
Using the matrix inversion lemma lfl4l . ( [27] ) can be written as 

(H^HQ* + J)- X H H 1EL - G H G(Q*G^G + I)" 1 = /jl — M . (30) 



Left multiplication by (H H HQ* + I) and right multiplication by (Q*G H G + I) of both sides of pO] ) 
yields 

H^H - G H G = [i (H^HQ* + I) (Q*G H G + I) - M (31) 
= // (G H GQ* + I) (Q*H H H + I) - M , (32) 
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where in obtaining (31 ) we have used the KKT condition (29 1, and Eq. (32 1 comes from the fact that (31 
is Hermitian. 

We are considering problem (25 1 for the case that H H H — G^G is strictly positive definite, i.e 



H^H y G H G, since this is the necessary condition for having a full rank optimal Q*. As we characterize 
the full rank Q*, the sufficient conditions are revealed as well. 

Remark 3. By following exactly the same steps as in the proof of [3, Proposition 5], one can easily 
show that for the case of H^H >- G H G, the optimization problem (25) is convex^] in Q. 



Thus for the case of interest, the KKT conditions ( [28] ), ( |29| ) and ( |3T| ) are necessary and sufficient 
conditions for the optimality of Q*. In other words, any Q >z that satisfies those conditions is an 



optimal solution for the problem (25 1. By the KKT condition (29), a full rank Q* follows that M = 0. 



Thus, (31 1 and (32 1 simplify to 

H H H - G^G = n (H H HQ* + I) (Q*G H G + I) 
= n (G^GQ* + I) (Q*H^H + I) 



(33) 
(34) 



Lemma 4. Let the diagonal matrix D and the unitary matrix 3>s respectively denote the eigenvalue and 
eigenvector matrices of (S^G^GS^ + 1), where we set S = (H H H - G H Gy 1 : 



8*G H G8* + I) = 3>s D 



H 



Then we have 



H H H = S 2 D 



G H G = S-s *g (D - I) 



Hi 



(35) 

(36) 
(37) 



Proof: Eq. (37 1 comes directly from (35). Please refer to Appendix A for details on obtaining (36) 



Using (36 1 and (37 1 in (33), after some simplification we have 



H (SH *s D *f S~5 Q* + I) (Q* SH *s (D - I) *f sH + I) 

(D *f S~5 Q* + Sf S>) (Q* S~* (D — I) + #f S~» 

H SH *s (D W + 1) *f S*s (W(D - I) + I) *f S~ i , (38) 



2 In fact, the optimization problem j^Jl is convex in Q when H H H >: G^G. 
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where in d38l) we denned 



From (38 1, we get 



I = /i (D W + I) *f S*s (W(D - I) + I) 



(39) 



(40) 



Let X >- and Y >- be two diagonal matrices, which will be defined soon. Left multiplication by 



X and right multiplication by Y of both sides of (40 1 gives 

XY = n (XD W + X) *f S*s (W(D - I) Y + Y) . 

We find X and Y by solving the following set of equations 

XY = I 

XD = (D - I)Y 

which results in 



X = Y" 1 = (I-D 1 ' 



Substituting ( |43| ) into ( |4T| ), we have 

I = /x(XDW + X) *fS*s (WDX + X" 1 ) 
= /x (XD WDX + X 2 D) X _1 D _1 *fS*gD _1 X _1 (XDWDX + D) 
= H (XDWDX + (D -I)) X- 1 D- 1 *fS* B D -1 X- 1 (XDWDX + D) 



where we have used (43 1 in obtaining ( 44 1. 

Using an approach similar to what we used to obtain (44 1 from (33 1, one can show thaj^] 



I = /x (XDWDX + D) X _1 D _1 *f S*iD _1 X _1 (XD W DX + (D - I)) 
Define K as 

K = XD W DX = (I - D 1 ) ' D W D (I - D _1 ) ' , 



where D and W are respectively given by (35 1 and (39 1. Moreover, let 



S = X _1 D _1 *f S* g D X X 1 = (I - D~ ) ~ 2 D" $f S*sD" (I - D _1 ) 



Clearly 1 45 i is also trivially obtained from ( 44 1 



(41) 



(42) 



(43) 



(44) 



(45) 



(46) 



(47) 
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Using ([46]) and ([47]» in (|44]> and (j45]>, we have 

I = /x (K + D) £ (K + (D-I)) 
= /x (K + (D-I)) £ (K + D) 



(48) 



We note from (48 1 that 



/i£ = (K + D)~ 1 (K + (D-I))" 1 = (K + (D-I))~ 1 (K + D) 



(49) 



This result implies that, for the optimal Q*, (K + (D — I)) -1 , (K + D) -1 and £ commute and have 
the same eigenvectors ifHl . 

Denote the eigenvalue decomposition of £ as 

s = unu ff . (50) 



Based on the argument made after (49), we have 

K + (D - I) = U Ai U H 



K + D = UA 2 U 



H 



(51) 



where one can easily confirm that A 2 = Ai + I. By replacing (50) and (51) in (48), and noting that 
U^U = UU^ = I, and that Ai, A 2 and ft are all diagonal, we have 



I = // UAi ft A 2 U H = n Ai ft A 2 = fi (Af + A x ) ft 



(52) 



Recall that the unknown parameters in (52) are the diagonal matrix Ai and the scalar // > 0. Let Xn 



and oji denote the ith diagonal element of Ai and ft, respectively. From (52), we can solve for Aji and 
obtain 



(53) 



where the Lagrange multiplier \i > is chosen to satisfy the power constraint Tr(Q*) = P t , as will be 
explained below. 

Theorem 2. The optimal full-rank input covariance matrix that attains the secrecy capacity for the average 
power constraint is given by 

Q* = S§* g (I - D -1 ) - * D 1 (UAiU H + I - D) D 1 (I - D -1 ) - * *f (54) 



. H^H - G^G >- 0, and 
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for the given P t , 



UAiU^ >- D - I 



(55) 



where Ai is given by (53 1 



Proof: The proof is obtained by obtaining K from (51 ), and substituting it back into (46) and (39 1 
via straightforward computations. ■ 
To fully characterize Q*, one must obtain the Lagrange multiplier \i > such that Tr(Q*) = P t . As 



(53 1 shows, Tr(Q*) is monotonically decreasing with /x: 

and 



lim Tr(Q* 



oo, 



lim Tr(Q* 

fj,— >oo 



-TVfSaSgD-^fSi) < . 



Thus for any transmit power P t , there exists a Lagrange multiplier \i > for which Tr(Q*) = P t . The 
appropriate value of \i can be easily found using, for example, the bisection method. 

Note that Theorem [2] also reveals the necessary and sufficient conditions for having a full rank optimal 
Q*. While for any transmit power P t one can find a Lagrange multiplier \i > for which Tr(Q*) = P t 



and UAiU H >- 0, to have a full rank Q*, (55 1 must be satisfied. Also recall from (74 ) that D - I y 



The flowchart in Fig. [T] summarizes the steps required to calculate the optimal full-rank Q* for the case 
of U H H - G H G y 0. 



Remark 4. By replacing Q* given by (54) into (25), the optimal input covariance matrix Q* given by 
Theorem [2] attains the secrecy capacity 

C sec {P t ) = log | A x | - log | Ai + I| + log |D| - log |D - I| 

= log|I-(A 1 +I)- 1 |+log|I + (D-I)- 1 | , (56) 



where Ai and D are diagonal matrices, respectively given by ([53]) and (35). Note that while both log 



terms in (56) return non-negative values, the first term depends on both the channels and the power P t , 
while the second term just depends on the channels (see Lemma [5]). 

In the following we show that for the case of H^H — G^G y 0, i.e. when at least one of the 
eigenvalues of H^H — G^G is zero, there exists an equivalent wiretap channel with the same secrecy 
capacity of the original wiretap channel. For this case, the secrecy capacity is characterized too. 

Theorem 3. For the MIMO Gaussian wiretap channel defined by direct and cross channels H and G 
respectively, with H^H — G^G y 0, there exists an equivalent wiretap channel with H eiJ and G eq such 
that Hf„H e9 - G"G eq y 0, and C sec (H, G, P t ) = C sec (H eq , G eq , P t ). 
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H, G, P t 

~~r~ 

Step 1 : Create S 



S = (H H H-G H G^ _1 



Step 2: Calculate #s and D from {35 

\ 



[*s, D] = eig(S2G ff GS2 +1) 



I 

Step 3: Calculate S from (47] 
1 



S = 



(I D 1 ) 2 D ^fS^D 1 (I D x ) 



Step 4: Calculate U and 17 from (50 



[U, O] = eig(S) 

J 

Step 5: Derive Ai from (53 
i 



A! = | f-i + (i + An- 1 ) 4 



Step 6: Calculate Q* according to (54 
i 



Q* = S2* § D _1 (I-D 



-l\-2 



UAiU fl + I-D) (I-D 



Step 7: Find 



M|Tr(Q*) = P t 
{ 

Step 8: Check the validity of Q* 
I 



UAiU w y D - I 



Fig. 1: Flowchart for obtaining the optimal full-rank Q*. 



Proof: Please refer to Appendix B for the proof and the characterization of the equivalent channels 
H eg and G eq . ■ 
Although at this time a proof is unavailable, we conjecture that the secrecy capacity of any general 
MIMO wiretap channel is given by an equivalent wiretap channel with H eq and G eq such that H^H^ — 
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G^ q G eq >- 0, and C sec (H, G, Pf) = C sec (H eq , G eq , Pt). If one cannot find such an equivalent wiretap 
channel, the secrecy capacity of the main channel is zero, i.e., C sec (H,G,P() = 0. Besides the case of 
H^H — G^G y 0, the case for which the optimal input covariance matrix is rank 1 is another case 
supporting this conjecture. For the rank 1 case, it is shown in [6] that the optimal input covariance matrix 
is given by Q* = P t \i\i H , where u is the normalized principal eigenvector corresponding to the largest 
eigenvalue Ai of the pencil (I + PtH^H , \ + P t G H G). For this case, if Ai > 1, the secrecy capacity is 
C sec (H, G, Pt) = log(Ai), and the equivalent wiretap channel is defined such that h^h e(? = u^H^Hu 
and g^g e(? = u^G^Gu. 

V. REMARKS REGARDING Q* 



This section discusses some interesting points regarding the optimal solution in (54). For the following 
observations, one can assume when required that both conditions for a full-rank Q* given in Theorem 
Q are satisfied. Let ji, i = 1, • • • ,rit, be the generalized eigenvalues of the pencil (H^H , G^G). 



Then from the definition [15], 7« = of, where <7j is the ith generalized singular value of (H , G). 



Lemma 5. The second term in the secrecy capacity expression (56 1 is only channel dependent and is 
equal to E^iM^?)- 



Proof: From (72i we have: 

= G H G Ss$s (D - I)" 1 D (57) 
= G fl GS^s(l + (D-ir 1 ) , 



where ( |57] ) comes from ( |73] ). Thus, from the definition 11141 . the generalized eigenvalue matrix of 
(H H H , G H G) is (I + (D - I)- 1 ), which completes the proof. ■ 
Note that the definition of singular values here is slightly different from what is given in [4]. Here <7j 



may be oo, while this is not the case in (4). More precisely, from (73 1 for the case of rank(G) = n e < rit, 
nt — n e diagonal elements of D are equal to one, as mentioned in Remark [5] The 0{ corresponding to 
di = 1 tends to oo. 

Lemma 6. In the high SNR scenario (P t — > oo) and for the case of rank(G) = n t , the asymptotic form 



of the exact secrecy capacity (56 1 is simply given by 



C sec = log |I + (D - I)- 1 ) = £ log(^ 2 ) . (58) 
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Proof: For P t — > oo, the Lagrange parameter satisfies fi — > 0, as mentioned after Theorem [2] 



Moreover, for the case of rank(G) = n t , the matrix 5] given by (47 1 will have finite-valued eigenvalues. 



Thus as n — > 0, the elements of the diagonal matrix Ai, given by (53 1, go to oo (\n — > oo). Consequently, 
the first log term in (|56|) disappears as Pt — > oo. ■ 



It is also interesting to consider the optimal solution in ( 54 1 for the case that the eavesdropper's channel 



is very weak, e.g. G = 0. For this specific case, the wiretap channel simplifies to a point-to-point MIMO 
Gaussian link, where the optimal input covariance matrix under the average power constraint is known 
to be &h (^1 — &h> an d is found via the standard water-filling solution, where unitary <&h and 

diagonal Ah are obtained from the eigenvalue decomposition H H H = <1?#A#<1>^. 



Lemma 7. For the case of G = 0, the optimal solution in ( 54 ) simplifies to the conventional water-filling 
solution, where 



Q* = f> H (^ 1 I-A H 1 )+ *f . 



(59) 



Proof: Using (72 1 and via simple calculations, we note that for any G, Eq. (54i can be rewritten as 
Q* = S2*g (I - D" 1 )^ D 1 UAiU H D 1 (I - D" 1 )^ *f S» - (H^H)" 1 . (60) 



From ( f73) , when G ->• then D ->• I. Next from ( [72] ), *s -> &h- Using these facts in (j47j), and after 
some straightforward calculations, we have £ -> (D - I^D^A^ 1 = ft, and U -»■ I in g). Using 



these in ( |60| ), we have 

Q* -»• * H (D - I)~ 1 D _1 A^ 1 Ai*f - (H^H)" 1 



<MD - I)" 1 A" 1 ( -I + (I + -1(D - I) A^ 



4 

— l 
A* 



*g - (H^H)- 1 



I - (H^H)" 1 = * H (m" 1 ! - A^)+ *g 



(61) 



where in obtaining ( 6 1 1 we used the fact that D — > I when G — > 0. 



VI. Numerical Results 

In the first example, we consider a MIMO wiretap channel with nt = n e 
matrices given by 

r 0.32 - 0.52? 0.83 + 1.15? 



2, nt = 3 and channel 



H 



0.51-0.26? 0.06-0.15? 
-0.11 + 0.8H 0.29 + 0.68? 



(62) 
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1.3 



Secrecy Capacity 

log (X^ 




/ y 

/ S 
/ / 
/ / 
// 
1 

1 




\ rank(Q*)=1 forP ( <2.8 
// 



10 20 30 40 50 

Transmit power P 



Fig. 2: Secrecy Capacity versus Pt for nt = n e = 2 and n r = 3. Solid curve represents secrecy capacity 
and dotted curve indicates the achievable secrecy rate using a rank-one input covariance matrix. 



G 



(63) 



0.03 - 0.70i -0.32 - 0.32i 
0.24 -0.1H 1.36 + 0.18i 
which satisfy H H H - G^G >- 0. Fig. g shows the secrecy capacity as a function of transmit power 
Pt. For comparison, the figure also depicts the achievable secrecy rate using the input covariance matrix 
Q = P t uu^, which results to R sec = log Ai, as shown in [4]-|6]. Note that in this example, the optimal 
Q* is not full-rank for P t < 2.8. 

In Fig. [3] we consider another example of the case of H H H -G H G >-0, here with 
nt = 4 and channel matrices given by 

0.89 + 0.54i -0.06 + 0.60i 0.48 - 1. Hi 

0.46 -0.44 + 0.80i -1.07 + 0.63i 

1.40-0.13? 0.17-0.82i 0.59 - OMi 

0.43 - 0.23i 0.03 + 1.35* 0.44 - 0.07i 



H 



(64) 
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/ 
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1 
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Transmit power P 



40 



50 



Fig. 3: Secrecy Capacity versus Pt for m = n e = 3 and n r = 4. Solid curve represents secrecy capacity 
and dotted curve indicates the achievable secrecy rate using a rank-one input covariance matrix. 



0.46 - 0.59i 0.24-0. OH -0.37 + 0.15i 

0.51-0.63? 0.58 + 0.5H 0.86 - 0.47i • (65) 
0.17-0.24i -0.83 + 0.51i 0.04 - 0.64i 
For this example, the optimal Q* is only full-rank for P t > 0.5. 

Finally in Fig. |4| we compare the standard point-to-point capacity without secrecy constraints with 



the secrecy capacity given by (56). In this example, P t = 20, direct channel H is given by (64i but 
the cross channel G is assumed to satisfy G^G = al, where a changes from to 1.95 (note that 
H^H — G^G y- only for a < 1.95). As predicted, the secrecy capacity achieved by the derived Q* 
in ( ]54| ) approaches the standard capacity as G — > 0. It is interesting to note that even for very small 
values of a, the difference between the standard capacity and secrecy capacity is considerable. 



VII. Conclusion 

In this paper, we considered the rank property of the optimal input covariance matrix under the average 
power constraint for a general MIMO Gaussian wiretap channel, where each node has an arbitrary number 
of antennas. We obtained necessary and sufficient constraints on the MIMO wiretap channel parameters 
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P=20 



■ Secrecy Capacity 

■ Capacity (G=0) 




Parameter a 



Fig. 4: Secrecy Capacity versus a for Pt = 20. Solid curve represents secrecy capacity and dashed curve 
indicates the point to point capacity. 



such that the optimal input covariance matrix is full-rank, and we presented a method for characterizing 
the resulting covariance matrix as well. 



Appendix A 
Proof of LemmaH] 

Define S = (H^H — G^G) -1 and apply the generalized eigenvalue decomposition on the pencil 
(SaH^HS^ + I , S^G^GS^ + I) to obtain the invertible generalized eigenvector matrix C and the 
diagonal generalized eigenvalue matrix Ag as 

C H [s^H^HSi + ij C = As (66) 

C H [s^G^GS^ +ll C = I . (67) 



By subtracting (67) from (66 1, we have 



C H C 



I . 



(68) 
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Note that from Lemma [5J we have A§ — I y 0. Thus, C must be of the form lfl4l 



(69) 



where $g is an unknown unitary matrix. In the following, as we continue the proof, $g is also charac- 
terized. 



By replacing (69 1 in (66) and (67), it is revealed that the unitary matrix 3>s represents the common set 
of eigenvectors for the matrices SaH^HSa + 1 and S2 G^GSa + 1, and thus both matrices commute. 
In particular, 



H 



StH^HSt +1 
S^G^GSs +1 



11 



* g = As(A, - I)" 1 = I + (Ag - I)" 1 
* g = (A, - I)" 1 . 



Defining D = (Ag — I) , from ( 70 )-(7 1 1 and via straightforward computation, we have 



H^H = S 2 D 



Hi 



G H G = S-2&- s (D-I) *fS"* 



(70) 
(71) 

(72) 
(73) 



which proves ( [36] ) and < \TT} . Substituting ( [68] ) in ( [67) , we also have 

C 



I = c 



StG^GSf +1 



or equivalently 



C^S^G^GS^C + As-I 



21- A= = C^S^G^GS^C. 



Remark 5. Since C^S^G^GS^C y 0, it results that 21 - A g y 0. Equivalently, by defining D 
(A g - I) -1 y 0, we have 



I-D 1 >-0 



(74) 



d - 1 y 0. 

Note that from (73 ), if G^G y 0, then D - 1 >- and vice versa. As we will observe in Theorem [2J to 
have a full-rank optimal input covariance matrix Q*, having a full-rank G^G is not required. While we 
assume throughout the paper and without loss of generality that the diagonal matrix D — I is invertible, 
for the case of rank deficient G^G one can follow the calculations in this paper assuming e > for 
zero-diagonal elements of D — I and letting e | at the end (see Lemma [7]). 
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Appendix B 
Proof of Theorem[3] 

We want to obtain QJ, the optimal input covariance matrix that attains the secrecy capacity for the 
case of H^H - G^G ^ 0. We note that rank(H H H - G^G) = m < n t . Hence, from Theorems [j] 
and |2j QjJ is rank-deficient. 



The right hand side of ( 25 ) can be rewritten as 



R(Q) = log |I + H H HQ| - log |I + G^GQ| 
= log |(I + H^HQ) (I + G^GQ)- 1 ! 
= log |I + (H H H - G H G) Q (I + G H GQ)~ 1 \ , 



(75) 



where Eq. (75 1 is obtained from the matrix inversion lemma [14] (I + A) 1 = I — A (I + A) . Let the 



eigenvalue decomposition of H^H — G^G to be denoted as 

Am 



U H H-G n G = * 



H. 







H 



(76) 



where A m >z is a diagonal matrix of size m x m. Using ( |76| ) in f75| ), we have 

A m 





R(Q) = log 



1 + * 



* H Q(I + G^GQ) 1 



log 



1 + 



A m 




(77) 



where in obtaining (77 1 we have used the facts that = q/^/ H = I and |I + AB| = |I + BA|. 

Define Q = ^ H Q^f and G = G^Sf, so that the optimization problem in (25 1 can be rewritten as 



max -R(Q) , 

QhQ, Tr(Q)=P t 



where 



R(Q) = log 



1 + 



A m Q 





Q(I + G H GQ) 



(78) 



From right-hand side of (|78|, we see that the optimal Q is of the form 

Q m Q 



Q 







(79) 
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where Q m ^ is of size m x m. Write G^G as 



G"G 



J 3 



(80) 



where Ji, J2 and J3 are of dimensions m x m, m x (n t — m) and (n t — m) x (n t — m), respectively. 



By substituting (79 1 and (80) into ( |78| ), we obtain 

fl(Q) = log 



^mQm 





I + JlQm 



log 



log 



1 + 



^mQm 



A m Q m (I + JiQ 




(1 + JiQ™)- 1 
-JfQ m (I + JiQ ■rr. 







I 



log 



I + A m Q m (I + JiQ m ) _1 

Q 1 

log |I + A m Q m (I + JiQ m ) _1 | 
log |I + (Am + Ji)Qj - log |I + JiQ 



R(Qr 



Using ( [81] ), the secrecy capacity is given by 

C S ec{Pt) 



max R(Q 

Q m hQ, Tr(Q m )=P t 



(81) 



(82) 



Problem (82i shows that the secrecy capacity of a wiretap channel with H^H — G^G ^ is equal 
to the secrecy capacity of an equivalent wiretap channel with 



A m + Ji 



G^G eg — Ji 



(83) 
(84) 



where A m and Ji are respectively given by (76 1 and (80 1. It should also be noted that for the equivalent 



channel, H^,H 



eq^eq 



g£g 



eq^eq 



A m y 0. Thus, the optimal Qj^ can be computed using Theorem [2j as 
long as the equivalent channel satisfies the second condition in Theorem [2] Finally, by substituting 
back into (1791) we obtain 



q: 



Qm Q 





H 



(85) 



which completes the proof. 
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